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Abstract 

The integral of the Wigner function over a subregion of the phase- 
space of a quantum system may be less than zero or greater than 
one. It is shown that for systems with one degree of freedom, the 
problem of determining the best possible upper and lower bounds on 
such an integral, over all possible states, reduces to the problem of 
finding the greatest and least eigenvalues of an hermitian operator 
corresponding to the subregion. The problem is solved exactly in the 
case of an arbitrary elliptical region. These bounds provide checks on 
experimentally measured quasiprobability distributions. 

The Wigner function has been much studied since its introduction 
not only in the context of quantum physics @, but also in signal processing 
0. For a quantum system in a pure state, the Wigner function carries the 
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same information as the wavefunction, up to an unimportant constant phase. 
In the case of a mixed state, it carries the same information as the density 
operator. 

An important property of the Wigner function, one of several properties 
which distinguish it from classical probability densities, is that its integral 
over a given subregion of phase-space may be negative or greater than one. 
Quasiprobability distributions which, according to quantum theory, corre- 
spond to Wigner functions, have been measured in recent experiments, for a 
variety of states of light and matter [|, |5], |7|) |9| ITOH, and negative values 
have indeed been observed. These experiments are probing the basic struc- 
ture and predictions of quantum mechanics in a new way, and the prospect 
of increasingly accurate experiments of this type adds greatly to the interest 
in, and importance of, the theory of the Wigner function. 

We consider the problem of determining the best possible bounds on the 
integral of the Wigner function over a given subregion of the phase-plane 
of any system with one degree of freedom. We show for any subregion of 
a rather general type that this problem reduces to the problem of finding 
the greatest and least eigenvalues of an hermitian Fredholm integral operator 
corresponding to that subregion. The problem is found to be exactly solvable 
for any elliptical or annular subregion, and the bounds are given explicitly in 
the case of the ellipse. These best possible bounds provide new information 
about the structure of the Wigner function, differing from known results such 
as best possible bounds (|?p on the values of the Wigner function itself, bounds 



on integrals of powers of the function [11], or bounds on various moments 
of the function [Q. In particular, the new bounds determine the degree 
to which the integral of any Wigner function over an elliptical subregion of 
the phase-plane can lie outside the interval [0, 1] which applies to classical 
densities. In principle, they therefore provide checks on experiments of the 
type to which we have referred, because they must be respected by any 
measured quasiprobability distribution consistent with quantum mechanics. 

As we shall show, appropriately chosen oscillator stationary states (or sin- 
gle frequency light modes) lead theoretically to the exact attainment of these 
upper and lower bounds. Such states are perhaps the easiest to establish 
experimentally, and it is just such states for which quasiprobability distri- 
butions have been measured in some of the experiments mentioned above 

In what follows, we consider systems with one degree of freedom, with 
a Cartesian coordinate q and its conjugate momentum p . Our results refer 
to the Wigner function considered at a particular instant, and are therefore 
independent of any particular dynamics. We work in dimensionless variables. 
Appropriate dimensional factors will appear in what follows if each coordinate 



2 



q there is replaced by q/L, each momentum p by Lp/h, each wavefunction ip 
by L 1 / 2 ^, each phase space area A by A/ft, and each Wigner function W by 
ftW, where L is a suitable constant with dimensions of a length. 

Given a normalized wavefunction if) corresponding to a pure state 
the Wigner function is defined as 

I r°° 

W^{q,p) = - i)*{q + x)i){q-x)e l2px dx. (1) 

7T J-oo 

Then [| 

Jw^dqdp=l, J[W i ,] 2 dqdp = — , (2) 

where T denotes the (g, p) phase-plane. 

For a mixed state, the density operator p is positive-definite and hermitian 
with unit trace, and typically can be resolved in the form 

P = ^Pi\A){A\, Pi>0, J2Pi = 1 > ( 3 ) 

i i 

where the states are orthonormal. The corresponding Wigner function 
has the form 

Wp = E» w *. ( 4 ) 

i 

where W 1 p i is the Wigner function corresponding to the pure state More 
generally, the sum in (Q) and could be replaced in part or whole by 
an integral, but this does not significantly affect the argument of the next 
paragraph. 

It follows from (^) and (^) that any bound on the Wigner function, or on 
its integral over a given subregion S of T, must hold for all possible mixed 
states if it holds for all possible pure states. For example, if 

/wW«.P)«U*»i for all (5) 
J s 

then for any p as in (|3]), 

/ W p (q,p)dqdp = ^Pi W^{q,p)dqdp 
J S i J s 

> 5>L = L. (6) 

i 

Since a pure state can be regarded as a limiting case of a mixed state, it then 
follows that best possible upper and lower bounds on the Wigner function or 
its integral, when considered over all pure states, must also be best possible 



3 



upper or lower bounds when considered over all mixed states, although a 
bound that is attainable over pure states may not in general be attainable 
over mixed states. Bearing this in mind, we restrict attention in what follows 
to pure states. 

Best possible bounds on the Wigner function itself are known 

--<W f (q,p)<- (7) 

7T 7T 

for all normalized ip, for all (q,p) G T. It is easily seen that W$ = ±1/tt at 
the point (q,p) if and only if 

^{q - x)e ipx = ±^j{q + x)e~ ipx for all x . (8) 

The problem of interest here is to find best possible bounds on the 
'quasiprobability functional' corresponding to the subregion S, defined as 

QslW^} = f W i ,(q,p)dqdp 
J s 

Xs(qiP)W^(q,p)dqdp, (9) 



where xs is the function with the value 1 on S, and the value on the 
complement of S. 

It follows at once from ([7]) and (|9|) that 

-— <Q s [wg<— , (io) 

7T 7T 

where As = J s dqdp is the of S. 

In order to obtain stronger bounds than flTOf) , recall that each real- valued 
function T(q,p) on T can be associated with an hermitian operator T such 
that 

(V,Tty) = / T(q,p)W^q,p)dqdp 1 (11) 



where (ipi,ip2) is usual scalar product of wavefunctions. Here T can always 
be written clS db Fredholm integral operator, 



{T^){x)= K T (x,y)^(y)dy } (12) 

J — oo 

with hermitian kernel given in terms of the real- valued function T(q,p) as 
K T (x,y) = ±-r T((x + y)/2,p)e ip ^dp. (13) 

Z7T J-oo 
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Consider now the case when T(q,p) = Xs(q,p)- Comparison of (Jffy and 
([TTD shows that 

QsiW^} = {1>,Ks4>), (14) 
(K s ip)(x) = / K s (x,y)^(y)d yi (15) 

J — oo 

tfsfoy) = -L f°° X s((x + y)/2,p)e L * x -vUp. (16) 

Z7T J-oo 

It follows at once from (fUj) that the extremal values of QsfW 7 ^] are de- 
termined by the eigenvalue problem Kgijj = Xip with Kg as in flI5|). In 
particular, 

inf Q 5 = A min , sup <3<j = X m ax , (17) 

where X m in and A max are the least and greatest eigenvalues of K$ respectively 
(or more generally, the infimum and supremum of the spectrum of Kg)- Thus 
the problem of interest now becomes the determination of X m in and A max . 

In order to proceed, suppose that the subregion S C V has the general 
form shown in Fig. 1. 



p = F 2 (q) ' 


p 

s 

3 Vj 






H 

p = Fx(q) 




c q 



Fig.l. A typical region S in the phase-plane. 

Here F\ and F 2 are real-valued functions defined for b < q < c, and satisfying 
F 1 (b) = F 2 (b), Fi(c) = F 2 (c), and F 2 (q) > F^q) for b < q < c. Each function 
need only be piecewise continuous, and b = — oo and/or c = oo is allowed. 
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For such a subregion, the characteristic function has the form 

Xs(q,p) = 
and the kernel (16) becomes 



1 b<q<c, F x {q) < p < F 2 {q) 
otherwise , 



1 rM^) . , > 
K s (x,y) = — / , e^dp 



2ir M(2±£) 



(19) 



2iri(x — y) 

for 2b < (x + y) < 2c, and otherwise. Note that the singularity at x — y is 
only apparent. Then flT5| ) becomes 

/ ^7 . mdy. (20) 

J2b-x 2m[x — y) 

More generally, the subregion S may consist of several nonintersecting 
parts Si, S2, ... of the same general type, even on overlapping q- intervals. It 
is easily seen that in such Ks = Ks 1 + Ks 2 + • • • However, in general 

[Ks 1 , Ks 2 ] 0) e ^ c -; so that the bounds associated with different subregions 
cannot be added. 

Note also that the extremal values of Qs and Qs' are the same if S is 
transformed into S' by a canonical transformation of T of the form 

q = otq + f3p + 7 , p' = lip + vq + p , (21) 

where a, (3, 7, /1, v and p are real constants satisfying ap — (5v = 1. In 
particular, the case of any circular or elliptical region of area na? can be 
reduced to the case of a circular disk of radius a, centred at the origin. 

In this case, the operator K$ (let K a denote it now) is given from (|20|) by 

(K a ^)(x) = 

2a-x sin[(x -y)Ja 2 -(x + ?/) 2 /4] 

7 s nv)dy, (22) 

-2a-x 7T(X — y) 

for — 00 < x < 00, and it is not hard to check that K a commutes with the 
simple harmonic oscillator Hamiltonian operator H defined by 

ff^) = -^+i 2 ^). (23) 
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This is explained by the fact that H generates transformations of the wave- 
function corresponding to rotations in the phase-plane, which leave the disk 
invariant. It follows that for every value of a the eigenfunctions of K a are 
the oscillator eigenfunctions 

il> n (x) = H n {x)e~ x2 l 2 , n = 0,l,... (24) 



where H n is the Hermite polynomial fll3| . 

According to (0), the eigenvalue X n (a) of K a corresponding to the eigen- 
function (0), must equal the total quasiprobability on the disk of radius 
a, as determined by the Wigner function W n (say) corresponding to that 
eigenfunction. Since it is known [Q, |T^] that 

W n (q,p) = (-l) n TX- l L n (2[p 2 + q 2 ])e-( p2+ ^ , (25) 
where L n is the Laguerre polynomial [JT3 1 , it follows that 

A n ( a ) = (-l) n / L n (2u) e~ u du . (26) 
Jo 



Thus A (a) = 1 - e~ a , Ai(a) = 1 - (1 + 2a 2 )e" a , A 2 (a) = 1 - (1 + 2a 4 )e" a 

4 
3 



A 3 (a) = 1 - (1 + 2a 2 - 2a 4 + \a G )e~ a \ etc. 



In summary: 

2a-x sin[(x -y)Ja 2 - (x + y) 2 /4] 

t s VnW dy 

-2a-x 7T{X — y) 

= \ n (a)^ n (x) , (27) 

with ip n as in (|24| ) and A„ as in (|26|). 

Fig. 2 shows the graphs of A n versus a for n — 0, 1, 2, 3, and also the 
graphs of A max and \ m in (bold lines). Note that A max (a) = Ao(a) = 1 — e~ a , 
whereas the graph of X min has the peculiar scalloped shape shown, because 
A m m( a ) = Ai(a) for < a < a 1; A min (a) = A 2 (a) for ai < a < a 2 , etc., 
where a± is the greatest value of a at which Ai(a) = A 2 (a), a 2 is the greatest 



value of a at which A 2 (a) = A3 (a), etc. Thus a\ — 1, a 2 = y (3 + \/3)/2, eic 
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Fig. 2. Left to right: graphs of A n for n — 0, 1, 2, 3, and also of A max , A 
(bold lines). 
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With the introduction of the appropriate dimensional factors, the result is 
that the integral of any pure-state or mixed-state Wigner function over any 
circular or elliptical region with area na 2 h in the phase-plane, lies in the in- 
terval [A m i n (a), A ma3; (a)], in contrast to the integral of any classical density, 
which lies in [0,1]. According to quantum mechanics, any quasiprobability 
distribution determined by quantum tomography (in particular) is described 
by a Wigner function [§, ||, [10| . For such a distribution, the quasiproba- 



bility on disks of various radii, centred on regions where the distribution is 
most negative, for example, could be estimated and checked for consistency 
against the theoretical bounds. Of course, experimental data are inevitably 
subject to noise for various reasons. While there are known techniques to 
allow for noise in the reconstruction of densities from more primitive data 
P, jlUH , this would obviously limit the power of the proposed check. A more 



subtle complication is that reconstruction algorithms may invoke quantum 
mechanical arguments || ||, jlGf , and any check would be satisfied trivially 



in a given case if these arguments forced a reconstructed density to satisy 
the theoretical bounds. Any given reconstruction algorithm would have to 
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be analysed carefully in this regard to ensure that a check was meaningful. 

If these difficulties could be overcome, might the proposed check be ele- 
vated to the level of a test of quantum mechanics itself? The assumptions 
underlying the theory of the Wigner function are very few: the linear vector 
space of states, the Born interpretation, and the conjugate relations between 
coordinates and momenta. However, quantum mechanics has been so well- 
tested at the energy scales of the present experiments that any violations, if 
indeed there are any, must surely be exceedingly small, and very probably 
beyond present capabilities of resolution amidst noise. 

Eigenvalue problems corresponding to other shapes such as squares and 
triangles are easily formulated, but do not seem to be exactly solvable. They 
could be tackled numerically. Exact results for disks can be extended to the 
case of an annular region (and more generally the case of several concentric 
annuli), because the operators K a commute for different a, and have common 
eigenfunctions. This may be particularly useful in checking distributions 
determined by the 'ring method' ||. 

These ideas can be extended to systems with more degrees of freedom, 
and to systems with spin. 

Thanks are due to J. A. Belward, R. Chakrabarti, G.A. Chandler, D. Elli- 
nas, W.P. Schleich and referees of a preliminary version for helpful comments. 
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